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Abstract We express connected Fermionic Green's functions in terms of two 
totally explicit tree formulas. The simplest and most symmetric formula, the 
Brydges-Kennedy formula is compatible with Gram's inequality. The second one, 
the rooted formula of [ARl], respects even better the antisymmetric structure of 
determinants, and allows the direct comparison of rows and columns which cor- 
respond to the mathematical implementation in Grassmann integrals of the Pauli 
principle. To illustrate the power of these formulas, we give a "three lines proof" 
that the radius of convergence of the Gross-Neveu theory with cutoff is indepen- 
dent of the number of colors, using either one or the other of these formulas. 



I) Introduction 

Perturbation theory for fermion systems is often said to converge, whether 
for boson systems it is said to diverge. But what does this mean exactly? Un- 
normalized Fermionic perturbation series with cutoffs are not only convergent but 
entire, whether Bosonic perturbation series with cutoffs have zero radius of conver- 
gence. But this is of little use since in actual computations, especially renormal- 
ization group computations, one needs to compute connected functions. It was not 
obvious until recently how to compute these connected functions in a way which 
preserves the algebraic cancellations of the Pauli principle to obtain convergent 
power series with accurate and easy estimates on the convergence radius, without 
using the heavy techniques of cluster and Mayer expansions. 
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The first solution to this problem which does not require cluster and 
Mayer expansions is in [FMRT] where some inductive way of computing such 
Fermionic connected functions is given and applied to the uniform convergence of 
two-dimensional many fermion systems in a slice. However in the corresponding 
computation the expansions steps were intertwinned with constructions of layers 
of some tree, and the final outcome was therefore not totally explicit yet. 

In [ARl] we developed two explicit tree formulas which are especially 
suited for cluster and Mayer expansions. Shortly thereafter we made the remark 
that the second formula leads to a completely explicit version of the inductive 
cluster expansion used in [FMRT] and [FKLT] . More recently we realized that the 
first formula is also compatible with Gram's inequality, and can therefore also be 
used in fermionic expansions. It is possible to develop a continuous renormaliza- 
tion group formalism for Fermions using this formula [DR], a problem which has 
attracted attention recently [S]. 

We also mention that an other (closely related) version of this kind of 
expansions has been introduced recently, using Gram's inequality [FKT] . It relies 
on a resolvent expansion rather than an explicit tree formula. 

II. Tree formulas and Grassmann integrals 

In this section we recall the notations and the "Forest Formulas" of [ARl] 
and apply them to Grassmann (fermionic) integrals. 

Let n > 1 be an integer, = {1, . . . , n}, Vn = G /„, i ^ j] 

(the set of unordered pairs in In). An element I of Vn will be called a iinif, a 
subset of Vn, a graph. A graph ^ = {/i, . . . , l^-} containing no loops, i.e. no subset 
{{^1) ^2}, {^2, ^3}, • • • , {^fc, ^1}} with /c > 3 elements, is called a forest. 

A Forest formula is a Taylor formula with integral remainder, which 

expands a quantity such as exp^^^^p^ ui^ to search for the explicit presence or 

absence of links ui . Taylor formulas with remainders in general are provided with a 
"stopping rule" and forests formulas stop at the level of connected sets. This means 
that two points which are joined by a link are treated as a single block. (More 
sophisticated formulas with higher stopping rules are useful for higher particle 
irreducibility analysis, or renormalization group computations but are no longer 
forests formulas in the strict sense [AR2]). Any such forest formula contains a 
"weakening factor" w for the links which remain undcrivcd. Many different such 
forests formulas exist, with different rules for w. Two of them were identified as 
the most natural ones in [ARl], corresponding to two different logics: in one of 
them, the Brydges-Kennedy formula, the forest grows in the most symmetric and 
random way, and in the other, the "rooted formula" it grows layers by layers from 
a preferred root. 

A forest is a union of disconnected trees T, the supports of which are 
disjoint subsets of In called the connected components or clusters of ^. 

The Brydges-Kennedy formula is then 
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Brydges-Kennedy Forest Formula 



u — forest 

(11.1) 

where the summation extends over all possible lengths r of ^, including r = 

hence the empty forest. To each link of ^ is attached a variable of integration wi; 
andw^.^y^ (w) — infjw;, / G L^{ij}^ where L^{ij} is the unique path in the forest 

^ connecting i to j. If no such path exists, by convention w^.^-^^ (w) = 0. 

The rooted formula is absolutely identical, but with a different rule for 
the weakening factor w, now called w^-^{w). It is a less symmetrical formula since 
we have to give a rule for choosing a root in each cluster. For each non empty 
subset or cluster C oi In, choose rc, for instance the least element in the natural 
ordering of In = {1, . . . ,n}, to be the root of all the trees with support C that 
appear in the following expansion. Now if i is in some tree T with support C we 
call the height of i the number of links in the unique path of the tree X that goes 
from i to the root rc- We denote it by l'^{i). The set of points i with a fixed 
height k is called the A;-th layer of the tree. The Rooted Forest Formula is then: 

Rooted Forest Formula 



\eVr, ^ 5={il,...,Jr} ^ = 1-^0 ^ ^V=\ ^ leVr^ 

u — forest 

(n.2) 

where the only difference with (II- 1) lies in the definition of w^.^y{w) which is 

different from the one of w^.^^^ (w) = 0. We define still w^.^{w) = if i and j 
are not connected by the ^. If i and j fall in the support C of the same tree 1 of 
^ then 

w^l^y{w) = if — I'^U)] > 2 (^z and j in distant layers) 

w^.^{w) — 1 if /^(z) = /^(i) (i and j in the same layer) 

wfy}(w) = w^ii>y if - 1 = = and {ii'] e T. (i and j 

in neighboring layers, i' is then unique and is called the ancestor of i in%). 

Proofs of these formulas are given in [ARl]. 

To any such algebraic formula, there is an associated "Taylor" tree for- 
mula, in which ui is interpreted not as a real number but as a differential operator 
ui = ^ [ARl] . Let us recall these associated Taylor formulas. Let S be the space 

of smooth functions from H^" to an arbitrary Banach space V. An element of 
H^" will be generally denoted by x = {xi)i^j, . The vector with all entries equal 
to 1 will be denoted by 1. Applied to an element H of S, the Taylor Rooted Forest 
formula takes the form: 
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The Taylor forest formulas 

E (nr'""')((nl:)^)w°"'(-)) ■ ("■3) 

5 u-forest ° ^ l€d ^ 

Here X^^ °' ^(w) the vector {xi)i^^^ o/M^" defined by xi = wf^^ ^(w), 
which is the value at which we evaluate the derivative of H. The symbol BK or R 
means that the formula is true using either the Brydges- Kennedy or the rooted 
weakening factor w. 

Again the proof of (11.3) is in [ARl] (but we recall that it is a rather 
trivial consequence of (II. 1-2) applied to ui = ^)- 

Any forest is a union of connected trees. Therefore any forest formula 
has an associated tree formula for its connected components. And therefore, at 
least formally, any forest formula solves the problem of computing normalized 
correlation functions. Indeed applying the forest formula to the functional integral 
for the unnormalized functions, the connected functions are simply given by the 
connected pieces of the forest formula, hence by the corresponding tree formula. It 
is in this sense that forest formulas exactly solve the well known snag that makes 
connected functions difficult to compute. This snag is that since typically there 
are many trees in a graph, one does not know "which one to choose" when one 
tries to compute connected functions in the (desirable) form of tree sums. Any 
forest formula gives a particular answer to that problem. It tells us exactly by 
how much our pondered "tree choice" has "weakened" the remaining loop lines: 
just by the weakening factor w. 




Figure 1 

To illustrate this by the simplest possible example, let us consider the 



second order connected graphs of the pressure for the (p theory in dimension 
with propagator 1, which simply counts multiplicities of connected graphs. There 
are two second order connected graphs, pictured in Fig 1. To count contractions 
correctly it is convenient to label all fields or half lines. The first graph Gi cor- 
responds to 2.Cl.Cl = 72 contractions, and has therefore weight 72. The second 
graph has weight 41—24. Now to select a tree in this case is to select the trivial 
trees T and there are 4^ = 16 such trees. If the same measure was used for the 
contractions of the three loop lines, we would find a weight 16.C|.C| =144 for 
Gi, and a weight 16.3! =96 for G2- This illustrates the overcounting problem of 
choosing a tree, and we see that this overcounting varies with the graphs. Now 
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we see how the weakening of the loop hnes corrects precisely this defect: with 
the weakening parameters w^^ or w^, the graph Gi is weakened by a factor 
Jg wdw = 1/2 and the graph G2 is weakened by a factor w^dw = 1/4 so that 
correct weights are recovered! The reader is invited to try higher order cases to 
see the distinction between w^^ and appear. 

It remains to see which formulas can be used in the constructive sense 
and for which theories. It is well known that the Brydges-Kennedy Formula 
[BK][AR1], when applied to interpolate a symmetric positive matrix, preserves 
positivity. Indeed wfj^^{w) completed as a symmetric matrix by the definition 

w^'^^(w) = IVi is a positive matrix, since it is a convex combination of block 
matrices with 1 everywhere [ARl]. Therefore if C is a positive matrix (like a 
bosonic covariance) the matrix wfj^^ {w)C{xi,Xj), being the Hadamard product 
of two positive matrices, remains positive for any values of the tree parameters w. 
This means that it is suited to the construction of bosonic models ( [B] [AR2] and 
references therein). But this positivity property also allows Gram's estimates for 
fermionic theories: 

Lemma 

Let A = {aaf3)a,/3 be a Gram matrix: a^p =< fa^dp > for some inner 
product < ., . >. Suppose each of the indices a and (3 is of the form {i, a) where 
the first index i, 1 < i < n, has the same range as the indices of the positive matrix 
{wfj^^ {w))ij, and a runs through some other index set S. 

Let B = {hctf3)a.,i3 be the matrix with entries bf^i^^^i^j^-j-^ = wfj^^ (w). < 
f{i,a)j 9{j,T) > (^''T'd let {bai3)aeA,i3eB be some square matrix extracted from B, then 
for any w we have the Gram inequality: 

\det{ha.p)a.^A,P^B\<\{\\fa\\X{\\9p\\ ■ (11.4) 



Proof: Indeed we can take the symmetric square root v of the positive matrix 
^S,BK gQ that wfj^^ = Xlfc=i '^ikVkj- Let us denote the components of the vectors 
/ and g, in an orthonormal basis for the scalar product < ., . > with q elements, 
by and g^ I < m < q. (Indeed even if the initial Hilbert space is infinite 

dimensional, the problem is obviously restricted to the finite dimensional subspace 
generated by the finite set of vectors / and g). We then define the tensorized 
vectors and G(j,r) with components F^.^^^ = VikfJ^^^^ and = Vjkg}]^^) 

where 1 < k < n and 1 < m < q. Now considering the tensor scalar product 
< ., . >T we have 

n q 

< F{i,CT),G{3,r) >T= Yl '"^kVjkf(^,a)9{j,T) = b{i,a){j,T) ■ (H.S) 

k=l m=l 

By Gram's inequality using the < ., . >t scalar product we get 

\detibaf3)aeA,f3eB)\ < n II/«IIt n \\9p\\t (11-6) 
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but - ELlEm=l(-^fi^))^ = ELlEm=l^i^fc(/(Ta))^ = 

'Wii Em=i(/(Ta))^ = WfMlfi since wu = 1 for any z, 1 < z < n. □ 
Let us now apply the Brydges-Kennedy formula to the computation of 
the connected functions of a Fermionic theory. The corresponding Grassmann 
integral is: 

^ J rf//c(^>)i^(V'a>a)e^^^<"'^«^ (II.7) 

where C is the covariance or propagator, P is a particular monomial (set of external 
fields) and S is some general action. We take as simplest example the massive 
Gross-Neveu model with cutoff, for which the action is local and quartic in a 
certain number of Fermionic fields. These fields are two-components because of 
the spin. In a finite box A the action is 

Sa = ^. [ dx{J2M^)M^)f (11.8) 

where a runs over some finite set of "color" indices. A is the coupling constant. 
The covariance C is massive and has an ultraviolet cutoff, hence in Fourier space 
it is for instance r]{p) /{p + m) where 77 is a cutoff function on large momenta. We 
only need to know that C is diagonal in color space C(x, a; j/, a') = for a 7^ a', 
and that it can be decomposed as 

C(x,y)= D(x,t)E(t,y)dt (II.9a) 



with 

\C{x,y)\ < K- — (II.9b) 

' ^ - {l + \x-y\)P ^ ' 

jD{x,t)\^dH <K ■ l\E{x,t)\'^dH < K (II.9c) 

for some constant K. 

This decomposition amounts roughly to defining square roots of the co- 
variance in momentum space and prove their square integrability. It is usually 
easy for any reasonable cutoff model. For instance if is a positive function, we 

can define D in Fourier space as (JIj^^2y/a and E as ^~(p2"^2)3/4^'' • 

It is the Fermionic covariance or propagator C =< 'ipa{x)'ipa{y) > which is 
interpolated with the forest formulas. We obtain, for instance for the pressure that 
the formal power series in the coupling constant is a sum over trees on {1, ....,n}, 
with 1, a distinguished vertex sitting at the origin to break translation invariance 
(similar formulas with external fields of course exist for the connected functions): 

Fermionic Tree Expansion 

p= lim -]-LogZ{A)= lim -]-( [ d//c(V', ^)e^^^'^°''^''^) 

A^oo |A| A^oo |A| J 
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oo A'' „1 

n=0 ai,...,a„,6i,. ..,&„ = ! T n ZeT ° 

da;i...da;n5(a;i = 0) JJ (C(a;j(/), a;j(/))(5i) x det{bai3)aeA,i3eB ■ (11.10) 

r/ie sum over the tti 's and hi 's are over the colors of the fields and antifields of 
the vertices obtained by expanding the interaction and of the form: 

i^ai {Xi)lpai {Xi)lpbi {Xi)i^bi (Xi) (H.H) 

with 1 < i < n. The sum over X is over all trees which connect together the n 
vertices at xi, . . . ,Xn- The sum over Q is over the compatible ways of realising 
the bonds I = {«, j} E % as contractions of a i/^ and '0 between the vertices i and 
j (compatible means that we do not contract twice the same field or antifield). A 
priori there are 8 choices of contraction for each I — {i-ij}, including the choice in 
direction for the arrow. e(T, O) is a sign we will express later. For any I e X, i{l) G 
{1, . . . , n} labels the vertex where the field, contracted by the procedure VL concerning 
the link I, was chosen. Likewise j{l) is the label for the vertex containing the 
contracted antifield. Si is 1 if the colors (among ai, . . . , a^, • • • , h^) of the field 
and antifield contracted by I are the same and else is 0. Finally the matrix {haj3)a,i3 
is defined in the following manner. 

The row indices a label the 2n fields produced by the n vertices, so that 
a = {i, a) with 1 < i < n and a takes two values 1 or 2 to indicate whether the 
field is the second or the fourth factor in (11.11) respectively. The column indices 
(3 label in the same way the 2n antifields, so that (3 — {j, r) with 1 < j < n and 
T = 1 or 2 according to whether the antifield is the first or the third factor in 
(11.11) respectively. The a's and j3's are ordered lexicographically. We denote by 
c{i,a) the color of the field labeled by (i, ex) that is ai, if a = 1, and hi if a = 2. 
We introduce the similar notation c{j, r) for the color of an antifield. Now 

b{i,T)(j,T) = w'^'^^{w)C{xi,Xj)6{c{i,a),c{j,T)) . (11.12) 

Finally each time a field (z, a) is contracted by 0, the corresponding row is deleted 
from the 2n x 2n matrix (hap). Likewise, for any contracted antifield the corre- 
sponding column is erased. A and B denote respectively the set of remaining rows 
and the set of remaining columns. The minor determinant featuring in formula 
(ILIO) is now det{boii3)a€A,/3&B which is (n+l) x (n+ 1). Indeed for each of the 
n — 1 links of %, a row and a column are erased. 

Although it is not important for the bounds, we indicate the rule for 
computing the sign e(T, O). Let (1,1) < cti < ••• < an-i ^ (^; 2) and 
(1,1) < (3i < ••• < (3n-i < {n,2) be the erased rows and columns respectively. 
Suppose that ar was contracted by T and Q together with P-n-{r), for 1 < r < n — 1. 
Clearly tt is a permutation o/{l,...,n — 1}. Then 

e{%,Q) = (-l)^(«i)+-+^(«"-i)+^(/'i)+-+^(/3"-i)e(7r) (11.13) 

where e(7r) is the signature of tt and for any a = (z, a) we introduced the notation 
i{a) = i and a{a) = a and similarly for any j3 = {j,T), we write j{j3) = j and 

t{P) = T. 
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Proof: The proof of the previous expansion although a bit tedious is straight- 
forward. First expand exp Sa{'4','4')- Each term will involve a Berezin integral 



/n 
diic{ip, V') n (^«* (^OV'ai {xi)'ipbi {xi)ipbi (xi)) (11.14) 
i=l 

which is a 2nx 2n determinant like the one of the matrix {hap)- Then in order to use 
(II. 3) we introduce for any pair {i, j} G a weakening factor xs^ij^ multiplying 
each of the eight entries involving both vertices i and j. The output of (II. 3) is 
a forst formula when applied to Z[K). However the amplitudes corresponding to 
each connected tree in the forest factorize (actually, to check that, one needs to be 
careful and compute e(T, O)). Taking the sum over trees instead of forests simply 
computes logZ'(A) instead of Z{K). Finally 5{xi = 0) accounts for the division by 
the volume |A| to get the pressure as an intensive quantity. □ 

Let us now use this fermionic tree formulas for proofs of convergence. 

Ill) Convergence of the tree formulas 

A typical constructive result for this Gross- Neveu model with cutoff is to 

prove: 
Theorem 

The pressure and the connected functions of the cut-off Gross-Neveu 
model are analytic in \ in a disk of radius R independent of N. 

Proof: The determinant det{hcii3)aeA,f3eB is diagonal in the colors of the in- 
volved fields and antifields. Each of the block determinants falls in the category 
described by the lemma, and thus the complete determinant is bounded by 
thanks to the decomposition (II. 9a). The spatial integrals are bounded using 
(II. 9b) for the propagators corresponding to the links in T, by K'^~^ . The sum 
over O is bounded by 8"^"^ (the number of compatible contractions is actually 
equal to 8 x 12"^'* x 10"^^ x 6"^^ where rip is the number of vertices of the tree T 
with coordination number p). The sum over colors is bounded by N'^'^^ , indeed 
once T and O are known, the circulation of color indices is determined. The at- 
tribution of color indices costs A^^ at the first vertex and by induction a factor 
A^ for each of the remaining vertices of the tree. Indeed climbing inductively into 
the tree layer by layer, at every vertex there is one color already fixed by the 
line joining the vertex to the root, hence one remaining color to fix, except for 
the root, for which two colors have to be fixed. The number of T's is bounded 
by n"^"^, by Cay ley's theorem. Finally the n-th term of the scries is bounded by 
(|A|'^/A^'^?i!)Ar"^+^i^'^~^8"^~-'^A""'+-'^n'^~^ and the radius of convergence is thus at 
least 1/86X2. n 
This is perhaps the shortest and most transparent proof of constructive 
theory yet! Remark in particular that it does not require any discretization of 
space, lattice of cubes, cluster or Mayer expansion, all features which are necessary 
for bosonic theories. This "three lines" treatment of the theory with cutoff can 
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presumably be extended into a "three pages" treatment of for instance the Gross- 
Neveu theory with renormahzation in two dimensions [DR] . Such a treatment with 
no discretization of space, no discrete sHcing of momenta, and continuous instead of 
discrete renormahzation group equations makes the constructive version of these 
fermionic theories almost identical to their perturbative version (and no longer 
more difficult). 

The theorem above is interesting not only for the analysis of the Gross- 
Neveu model but also for that of the two-dimensional interacting Fermions consid- 
ered in [FMRT] or [FKLT]. In this latter case, the "colors" correspond to angular 
sectors on the Fermi sphere and the factor 1/N in the coupling is provided by 
power counting. 

We now give a second and longer proof of the theorem using the rooted 
formula. This is worth the trouble since the weakening factor in the rooted formula 
completely factorizes out of the determinant. This second formula may therefore 
be useful in problems for which Gram's inequality is not applicable and the method 
of "comparison of rows and columns" of [IM] or [FMRT] has to be used. 

We return to formula (11.10) but use the rooted weakening factor. To 
bound the determinant, we need to introduce simply one further subtle modifica- 
tion. We need a further sum which decides for each field or antifield whether it 
contracts or not to the level in the tree immediately below it. 

The only change in formula (11.10) is that the loop determinant has en- 
tries 

W = "^»x5j(/3)(w)C(a;i(c,),a;j(^))5(c(Q!), c{(3)) (III.l) 
with q: e ^, /? e -B. Any of these entries we write as: 

with the convention that 

h%^h^p5{l'^{i{a)),l'^{3m) (ni.3a) 

= ba^Sil'^iiia)) + 1, l^UiP))) (111.3b) 

ft = ba^^iiicx)) - 1, l^Um) (ni.Bc) 

and blf^''^^^ = whenever {u{a), v{p)) = % ft), (ft, \\), (b, \\), % b), (b, b) or (U, U). 

Then we expand the determinant det{bai3)aeA,/3eB by multilinearity with 
respect to rows, having for any a E A, a, sum over u{a) G {ft, t], b}. Then we 
do the same thing for columns and we get a sum over v{P) e {tt^tljb}, for any 
P E B. We then obtain an expansion into 32"+2 determinants, of the form 

detibl^^^''^^\eA,f3eB. Let / = A,u{a) = b} and / = e B,v{(3) = b} 

be the number of flat fields and antifields. The number of sharp fields and antifields 
s, s, and of natural ones n, n are defined similarly. For one of the determinants 
to be nonzero we must have n = n, f = s and s — f. We can also group together 
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fields and antifields according to their level in the tree T, so that we get a block 
diagonal determinant that can be factorized as: 

kmax 

k=0 

X li ^^*(C)a.4,^.Bt,, ff ^^^(C)a€4,/3.B« , (111-4) 
k^O k=l 

where kmax is the maximal occupied level of the tree T. A\ is the number of cc's 
in A with l'^{i{a)) — k and u{a) — tj. A^, A\., b\, bI and B^. are defined in the 
same manner. We must have of course = ^(S^), ^(A^) = and 

i^{A\.) = H^{b\_^) for any /c, in order for the full determinant to be nonzero. Now 
by definition of w^^- (w) we have 

^e*(&a/3)c,eA5./3eBS. =det(c(xi(c,),Xj(^))5(c(Q;),c(/3)) ) . (III.5) 



Likewise 



= n ^j(/3)i'(^) <^(^i(a)'^j(/3))^(c(«)>c(/3)) ) (III.6) 



where denotes the ancestor of in the rooted tree T. Finally 

= Yl Wi^a)i'(a) X d(it\C{Xi(a),Xj(^P))S{c{a),c{P))] . (III.7) 

a€Ai ^ ^^eAlPeBl_^ 

k 



We see that the w's factor out and are bounded by 1. The remaining determinants 
again can be factorized according to the colors c{{a), and c{(3). Each determinant 
obtained that way is of the form det{C{xi(^a):Xj{f3)))ap- 

At this stage we can complete the proof of the Theorem either using 
Gram's inequality or the method of comparisons of rows and columns to bound 
these determinants. 

Gram's inequality is the fastest road. We use (11.9) to bound every p by 
p determinant by K^. And the complete loop determinant is then bounded by 

\detMaeA,/3€B\ < 32-+2K"+i . (III.8) 

The end of the proof is as in the case of the Brydges-Kennedy interpolation. □ 
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It is also clear on the form (III. 3) that the method of comparisons of rows 
and columns of [IM], [FMRT] also works here. This method roughly corresponds 
to Taylor expanding around a middle point further and further when fields or 
antifields accumulate in any given cube of unit size of a lattice covering IR'^. 
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